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Abstract. We study one of the central open questions in one-dimensional renormaliza- 
tion theory - the conjectural universality of golden-mean Siegel disks. We present an 
approach to the problem based on cylinder renormalization proposed by the second au- 
thor. Numerical implementation of this approach relies on the Constructive Measurable 
Riemann Mapping Theorem proved by the first author. Our numerical study yields a 
convincing evidence to support the Hyperbolicity Conjecture in this setting. 



1. Introduction 

One of the central examples of universality on one-dimensional dynamics is provided by 
Siegel disks of quadratic polynomials. Let us consider, for instance, the mapping 

Pg(z) = + e^'^'^z, where 9= {y/5 + l)/2 

is the golden mean. By a classical result of Siegel, the dynamics of Pg is linearizable 
near the origin. The Siegel disk of Pe, which we will further denote Ag is the maximal 
neighborhood of zero in which a conformal change of coordinates reduces Pg to the form 
w I— i> e^'^^^w. By the results of Douady, Ghys, Herman, and Shishikura, the topological 
disk extends up to the only critical point of Pg and is bounded by a Jordan curve. 

It has been observed numerically (cf. the work of Manton and Nauenberg jMNj ). that the 
boundary of is asymptotically self-similar near the critical point. Moreover, the scaling 
factor is universal in a large class of analytic mappings with a golden-mean Siegel disk. In 
1983 Widom |Wij defined a renormalization procedure for Pg which "blows up" a part of 
the invariant curve dAg near the critical point, and conjectured that the renormalizations 
of Pg converge to a fixed point. In addition, he conjectured that in a suitable functional 
space this fixed point is hyperbolic with one-dimensional unstable direction. 

In 1986 MacKay and Persival |MPj extended the conjecture to other rotation numbers, 
postulating the existence of a hyperbolic renormalization horseshoe corresponding to Siegel 
disks of analytic maps, analogous to the Lanford's horseshoe for critical circle maps [Lanll 
|Ean2j. 



In 1994 Stirnemann |Stirj gave a computer- assisted proof of the existence of a renor- 
malization fixed point with a golden-mean Sieg el disk. In 1998, McMuUen |McMj proved 
the asymptotic self-similarity of golden-mean Siegel disks in the quadratic family. He con- 
structed a version of renormalization based on holomorphic commuting pairs of de Faria 
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jdFH ldF2j , and showed that the renormahzations of a quadratic polynomial with a golden 
Siegel disk near the critical point converge to a fixed point geometrically fast. More gener- 
ally, he constructed a renormalization horseshoe for bounded type rotation numbers, and 
used renormalization to show that the Hausdorff dimension of the corresponding quadratic 
Julia sets is strictly less than two. 

Having thus attracted much attention, the hyperbolicity part of the conjecture of Widom 
for golden-mean Siegel disks is still open. 

In |Yalj the second author has introduced a new renormalization transformation TZcyi, 
which he called the cylinder renormalization, and used it to prove the Lanford's Hyperbolic- 
ity Conjecture for critical circle maps. The main advantage of TZcyi over the renormalization 
scheme based on commuting pairs is that this operator is analytic in a Banach manifold of 
analytic maps of a subdomain of C/Z. It is thus a natural setting to study the hyperbolic 
properties of a fixed point. In the present paper we study the fixed point of the cylinder 
renormalization numerically, and empirically confirm the hyperbolicity conjecture, as well 
as study the dynamical properties of the fixed point. The main numerical challenge in 
working with cylinder renormalization is a change of coordinate involved in its definition. 
It is defined implicitly, and uniformizes a dynamically defined fundamental domain to the 
straight cylinder C/Z. To handle it, we use the Constructive Measurable Riemann Map- 
ping Theorem developed for numerically solving the Beltrami partial differential equation 
by the first author in |GaiL IGKj . 

2. Definition and main properties of the cylinder renormalization of 

Siegel disks. 

Some functional spaces. For a topological disk W G C containing and 1 we will 
denote Aw the Banach space of bounded analytic functions in W equipped with the sup 
norm. Let us denote Cw the Banach subspace of Aw consisting of analytic mappings 
h:W such that h{0) = and h'{l) = 0. 

In the case when the domain W is the disk Dp of radius p > 1 centered at the origin, we 
will denote Ad^ = Ap and = Cp. 

For each p > 1 we will also consider the collection of analytic functions f{z) defined 
on some neighborhood of the origin with /(O) = 0, equipped with the weighted h norm on 
the coefficients of the Maclaurin's series: 



(2.1) ll/llp = E^4r^^'^- 

n=0 

We will further denote the subset of consisting of maps / with the normalizing 
condition /'(I) = 0. 

The proof of the following elementary statement is left to the reader: 

Lemma 2.1. 

1) LetfeLl, then sup„ Jf{z)\ < \\f\l; 



CYLINDER RENORMALIZATION OF SIEGEL DISKS. 



3 




Figure 1. Schematics of cylinder renormalization. 



2) Let f e Ap, and p' > p, then \\f\\p < ^ supp^, \f{z)\. 
As an immediate consequence, we have: 
Corollary 2.2. is a Banach space. 

Cylinder renormalization operator. The cylinder renormalization operator is defined 
as follows. Let / G Cw Suppose that for n G N there exists a simple arc / which connects 
a fixed point a of /" to 0, and has the property that /"(/) is again a simple arc whose 
only intersection with / is at the two endpoints. Let Cf be the topological disk in C \ {0} 
bounded by / and /""(/). We say that Cf is a fundamental crescent if the iterate /^"|cy^ 
mapping /"(/) to / is defined and univalent, and the quotient of Cf U f~^{Cf) \ {0, a] by 
the iterate /" is conformally isomorphic to C/Z. Let us denote Rf the first return map 
of Cj, and let us denote z the critical point of this map (corresponding to the orbit of 
0). Let g be the map Rf becomes under the above isomorphism, mapping z to 0, and 
h = e o g o e~^, where e{z) = exp [—2niz]. We say that / is cylinder renormalizable with 
period n, if /i G Cy for some V, and call h a cylinder renormalization of f (see Figure 1). 
We summarize below the basic properties of cylinder renormalization proven in |Yalj : 

Proposition 2.3. Suppose f G Cw is cylinder renormalizable, and its renormalization hf 
is contained in Cy- Denote Cf the fundamental crescent corresponding to the renormal- 
ization. Then the following holds. 

• Every other fundamental crescent Cf with the same endpoints as Cf, and such that 
Cf U C f is a topological disk, produces the same renormalized map hf. 

• There exists an open neighborhood U{f) C Cw such that every map g G U{f) is 
cylinder renormalizable, with a fundamental crescent Cg which can be chosen to 
move continuously with g. 

• Moreover, the dependence g ^ hg of the cylinder renormalization on the map g is 
an analytic mapping Cw Cy. 

We now want to discuss the dynamical properties of the cylinder renormalization of 
maps with Siegel disks derived in jYa3j . To simplify the exposition let us specialize to the 
case when the rotation number of the Siegel disk is the golden mean e = (V5 + l)/2. The 
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golden mean is represented by an infinite continued fraction 

9=1 + 5^=1 + [1,1,1,...]. 

1 + — 

As is customary, we will denote Pn/Qn the n-th convergent 

Pn/qn = [1,1,1,..., 1]. 

^ V ' 

n 

Theorem 2.4 ( |Ya3j ) . There exists a space Cu and an analytic mapping f G Cu which 
has a Siegel disk Aq with rotation number 9 whose boundary is a quasicircle passing through 
the critical point 1, such that the following holds: 

(I) There exists a branch of cylinder renormalization with period qt, k E N, which we 
denote TZcyi such that 

Tlcyif = /; 

(II) the quadratic polynomial Pe{z) = e^'^'^^z + z"^ is infinitely cylinder renormalizable, 
and 

T^^cyiPe /, 

at a uniform geometric rate; 
(111) the cylinder renormalization TZcyi is an analytic and compact operator mapping a 
neighborhood of the fixed point f in Cu to Cu- Its linearization C at f is a compact 
operator, with at least one eigenvalue with the absolute value greater than one. 

A central open questions in the study of TZcyi is the following: 

Conjecture 2.5. Except for the one unstable eigenvalue, the rest of the spectrum of C is 
compactly contained in the unit disk. 

Our numerical study of TZcyi will begin with empirically establishing the convergence to 
/. We will then make explicit the choice of the neighborhood U in the above Theorem. 
Experimental evidence suggests that it can be taken as a round disk Dp for some particular 
value of p. Having numerically established this, we will then proceed to experimentally 
verify the Conjecture. 

3. Construction of the conformal isomorphism to the cylinder 

The principal difficulty in numerical, as well as analytical, study of cylinder renormal- 
ization is the non-explicit nature of the conformal isomorphism 

^■.CfU f-{Cf) \ {0, a} ^ C* 

of a fundamental crescent, which is a part of the definition of TZcyi- An analytic approach 
to this construction based on the Measurable Riemann Mapping Theorem was presented 
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by the first author in \Ga.i\ . It has its roots in the complex-dynamical folklore; similar 
arguments are found, for instance, in the work of Lyubich |Lyu| and Shishikura |Shij . 

In |Gaij , the first author demonstrates how this approach can be implemented construc- 
tively, with rigorous error bounds. We will give a brief outline here. 

Uniformization of the cylinder using the Measurable Riemann Mapping The- 
orem. We will start with a description of our choice of a fundamental crescent with 

period qn for a map / G Cu sufficiently close to /. 

To construct the boundary curve /„ of C" consider first the union /„ of two parabolas 
x + i{Ax'^ + Bx) and {Cy'^ + Dy + E)+iy: the first passing through points and fi"+^~^i"(^l^^ 
the second — through yj^+a+^^^i) and a repelling fixed point a^^. All parameters in 
these two parabolas are defined uniquely after one specifies their common tangent line at 
fi"+2+i"(^l^ (see equations ()3.4j) below). While somewhat arbitrary, this choice has the 
virtue of possessing a simple analytic form. It can be shown rigorously (see jGaij ) that 
by modifying in sufficiently small neighborhoods of the endpoints (small enough not to 
influence our numerical experiments) we obtain a curve In which together with f~'^"{ln) 
bounds a fundamental crescent for TZcyi- 

Now consider the following conformal change of coordinates for z G CJ: 

The normalizing constant (3 will be chosen so that 

r-i(/^"+Hl)) = 0, 
while a real positive a will be specified by the condition 

\r~\p-+^+<i-{l))\ = 1. 

The choice of of this coordinate is motivated by the fact that maps the interior of the 
fundamental crescent conformally onto the interior of an infinite vertical closed strip S, 
whose width is comparable to one (cf. Figure 2). Next, similarly to PHij, define a function 

Qn-.K = {u + iv eC : < Rew < 1} — yS 

by setting 

gn{u + iv) = (1 - M)r-i(/-^"(7„(t;))) + ur-\-fniv)), 
where 7„ is a parametrization 

7„ : M ^ 

which we will specify below. 

Let (Jo be the standard conformal structure on C, and let a = g^cro be its pull-back on 
U. Extend this conformal structure to C through 

a = {T^Ya on T~^{U), where T{w) = w + 1, for all G N. 



6 DENIS GAIDASHEV, MICHAEL YAMPOLSKY 

Assuming the mapping gn is quasiconformal, the dilatation of a is bounded in the plane. 
By Measurable Riemann Mapping Theorem (see e.g. |ABj ) there exists a unique quasicon- 
formal mapping ^ : C i-^ C such that g*crQ = a, normalized so that ^(0) = and ^(1) = 1. 
Notice that g oT o g^^ preserves the standard conformal structure: 

and therefore it is a conformal automorphism of C. Liouville's Theorem implies that this 
mapping is affine. By construction, it does not have any fixed points in C, and hence, is a 
translation. Finally, g oT o ^"^(0) = 1, and thus 

g oT o g^^ = T. 

By the definition of (jn, 

9n'°r-'°r"°r = Tog—^ 

on the image of f~''"{ln) by r~^. Set (p = g o and $ = 00 r^^. Clearly, $modZ is a 
desired conformal isomorphism 



CjUf~HC])\{0,a,J—.C/Z. 

Again, set e{z) = e"^'^*^, and g = eo g o e~^. Since 

9zie{w)) ^ e{w) g^{w) 
9zie{w)) e{w) 9wiw)' 
the 1-periodic function ^ is a solution of the Beltrami equation 

gw = P'Ow, /i = {gn)u,/ {gn)w 

whenever g is a solution of 

(3.2) g^ = ng.,, n{z) = {z/ z)p,{e'\z)). 
Thus, we have reduced the problem of finding 

$ = eo<|> = (y(oeo g~^ o r~"^ 
to that of finding the properly normalized solution of the Beltrami equation 

(3.3) 9, = fig., M^) = ^7m-T73T 

on the punctured plane C*. 

It remains to describe the choice of the parametrization of in the definition of cjn- It 
is convenient for us to parametrize /„ using the radial coordinate in C For n = 1 and 
/ G Cui sufficiently close the empirical fixed point of the cylinder renormalization with 
period 1 we use the following parametrization: 

r {x{r),Ax{rf + Bx{r)), r <f, 
'^^^ \{Cy{ry + Dy{r) + E,y{r)), r > r, 
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where 



^[r) 1/4(1)1 M-^' 

?/(r) = Im/ (1) ^ — hlmai- 



Tir) 



\ai-m\ + \m 



and f is defined through T(f) = |/^(1)|. Constants A, B, C, D and E are fixed by the 
conditions 0, /^(l),ai G li, together with the requirement that the slope of the common 
tangent fine to both parabolas at the point /^(l) is equal to 1.1. 

This particular choice of the parametrization is motivated by the speed of convergence 
of the iterative scheme in the Measurable Riemann Mapping Theorem. 

We define the following function on C*: 

(3.5) 9^{r, 0) = (r/(-0) + ^) T-\r\X,{r))) + (l - .^(-0) - r-\X,{r)), 

where — vr < < vr, and rj is the Heaviside step function (we have adopted the convention 
7/(0) = 1). Then, according to ()3.3j) . the Beltrami differential fi is given by the following 
expression: 

(3.6) fi(re"^) = ^2^4. ^^r9l{rA) + id^gijr,^) 

rdrgiir, (p) - id^gi{r, ^) 

on C\(— oo, 0] . This is the expression that we have used to compute the Beltrami differential 
in our numerical studies. 

A constructive Measurable Riemann Mapping Theorem (MRMT). To solve the 
Beltrami equation numerically, we use the constructive MRMT proved by the first author 
in |Gaij . Before formulating it, we need to recall two integral operators used in the classical 
approach to the proof of MRMT (see |ABj ) . 
The first of them is Hilbert Transform: 

(3.7) T[h]{z) = ^\im J I jP^d^^d^, 
the second is Cauchy Transform 

(3.8) P[h]{z) = ^ [ [ j^.d^Adt 

Hilbert Transform is a well-defined bounded operator on Lp(C) for all 2 < p < oo. For 
every such p there exists a constant Cp such that the following holds (cf. |CZp : 

||r[^]|L < Cp \\h\\ for any h G ivp(C), and Cp — 1 as p — > 2. 
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We are now ready to state the constructive Measurable Riemann Mapping Theorem of 

Theorem 3.1. Let fi G Loo(C) and an integer p > 2 be such that < K < 1 and 

Kcp < 1, where 

Cp = cot^(7r/2p). 

Assume that = z/ + 77 + 7, where v and rj are compactly supported in D/j, and '~f{z) is 
supported in <C \ D/j. Furthermore, let rj be in Lp(D^) and Hr^H^ < 5 for some sufficiently 
small 5. Also, let h* G Lp{C) and e be such that Bp{h*,e), the ball of radius e around h* 
in Lp(C), contains BpiTlu^h* + l)],Cpe'), with 

e' = 5 essupn^l/i* + 1| + Ke. 

Then the solution g'^ of the Beltrami equation g^ = fJ^g^ admits the following bound: 

(3.9) \g^iz)-g:iz)\<Fie',R;z,g:iz),p,K,Cp) 

where g^{z) = P[u(h* + l)](;z) + z, and F{e', R) = 0(e', R~'^^p) is an explicit function of its 
arguments. 

Given the theorem, the algorithm for producing an approximate solution of the Beltrami 
equation is as follows. Given a /i as in the condition of the theorem, we first iterate 

(3.10) h^T[u{h + l)] 

to find a numerical approximation h* to the solution of the equation T[fi{h* + 1)] = h*. 
After that, we compute an approximate solution as 

(3.11) g:{z) = P[u{h: + l)]{z) + z. 

One can obtain rigorous computer-assisted bounds on such solution using Theorem 13.11 
Such bounds have been indeed implemented in |Gaij for a particular case of the golden 
mean quadratic polynomial. However, in the present numerical work we will not require 
such estimates. 

In the Appendix, we will discuss several numerical algorithms for the two integral trans- 
forms appearing in this scheme. 

4. Empirical convergence to a fixed point 

An appropriate choice of the domains of analyticity for the renormalized functions is cen- 
tral to a successful numerical implementation of cylinder renormalization. Our numerical 
approximation to the renormalization fixed point is a finite-degree truncation of a func- 
tion analytic in D3 (see Section |S1 for a detailed explanation of this choice of the domain). 
However, for the purposes of obtaining bounds on higher-order terms, we will consider a 
smaller analyticity domain, a disk of radius p = 2.266. Thus the cylinder renormalization 
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Figure 2. Schematics of renormalization. The contours (7(7) 1— > g{r) are 
used to find a polynomial approximation of TZcyif through the Cauchy inte- 
gral formula. 

will be a priori an analytic operator in a neighborhood of the fixed point in the Banach 
space Lp with p = 2.266. ^ 

Given an / G L^, a numerical approximation to its cylinder renormalization of order n 
is built as follows. As the first step, we construct a fundamental crescent as described 
above, and find the normalized solution g of the Beltrami equation 

gz{z) = fi{z)gz{z) with p as in ()3.6p 

as described in Section El and the Appendix. Next, we choose a contour 7 in the domain 
of g, and map this contour into the fundamental crescent hj r o g^. 

7 = r 0^1(7). 

Applying the first return map to the points of this contour, we obtain T = Rfij), and 
find the images, g{'y) and g{T) (where T = (jff ^(r^^(r))). The coefficients in a finite order 
polynomial approximation to 

Ticylf = 90 o o Rj o T o g^ o g-^ 

^We have implemented the procedure for the cyhnder renormahzation described in Section 13 and a par- 
ticular method of solving the Beltrami equation (see Appendix) as a set of routines in in the programming 
language Ada 95 (cf j^DAl for the language standard). We have parallelized our programs and compiled 
them with the public version 3.15p of the GNAT compiler ADA2'' . The programs ( |Prog| ) have been run 
on the computational cluster of 92 2.2 GHz AMD Opteron processors located at the University of Texas 
at Austin. 
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are then found via the Cauchy integral formula using these two contours ^'(7) and g(T) 
(see Figure 121) ■ 

As seen in Theorem 12.41 the sequence of the cylinder renormalizations of the quadratic 
polynomial Pe converges to a fixed point / = TZcyif- We have used this fact to compute 
an approximate renormalization fixed point fa as the cylinder renormalization IZ^yiPe of 
order k = 11. 

Further, we improved this approximation by iterating 

(4.1) fa^^sOTZ^yifa 

where is the projection on the candidate stable manifold of / 

W = {feLl: /'(O) = e'-''}, 

defined by setting 

In this way we have obtained a polynomial fa of degree 17, and have estimated that, 
not taking into account the errors in the solution of the Beltrami equation and due to 
round-off, 

(4.2) \\n,yifa - faWp < 1.88 X 10-3 ^ 0.89 x 10-^||/a||p. 

Moreover, the iteration 1)4.11) does not lead to a significant variation in the computed values 
for the coefficients of fa, which indicates that the original approximation is indeed quite 
accurate. The largest change is in the highest coefficient, which differs by 0.4% for fa and 
its renormalization. Of course, this represents a negligible correction to the absolute value 
of the coefficient itself. 

The approximate expression for fa is as follows (all numbers truncated to show six 
significant digits): 



/a(x) = 



X 














( 8.00882 X 10" 


-1 


+ Z4.07682 X 10" 


-i\ 


+ 


x^ 


(-4.12708 X 10' 


-1 


+ i2.97670 X 10" 


-2\ 


+ 


4 

X 


( 1.02033 X 10" 


-1 


- Z9.83702 X 10" 


-2\ 


+ 


x^ 


( 2.61573 X 10" 


-5 


+ i4.13871 X 10" 


-2\ 


+ 


x^ 


(-8.42868 X 10" 


-3 


- Z6.96474 X 10" 


-3\ 


+ 


7 

X 


( 2.60095 X 10" 


-3 


- i6.58544 X 10" 


-4\ 


+ 


x^ 


(-2.01382 X 10" 


-4 


+ i5.95113 X 10" 


-4\ 


+ 


Q 

X 


(-9.40057 X 10" 


-5 


-il. 11237 X 10" 


-4\ 


+ 




( 3.21762 X 10" 


-5 


- i4.40144 X 10" 


-6\ 


+ 
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5. Domain of Analyticity of the Renormalization Fixed Point 
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Compactness of TZcyi- We have verified experimentally the compactness property of TZcyi 
stated in Theorem 12.41 More precisely, we observe the following empirical fact: 



Set p = 2.266 and p' = 3. Then we can take U = 'Bp in Theorem \2.4\ More specifically, 
the fixed point f is a well-defined analytic mapping in Cp, and moreover, if we denote 

9 = fUp, then n^yig e Cp/. 

To verify the claim numerically, we have used the approximation fa obtained in the 
previous section. To estimate p', we have chosen a curve 7 in the fundamental crescent, 
such that $^^(7) is a simple closed loop that encircles D3 (Figure 3). We then verify that 
the orbit under the return map of the component of the set Cf \7, such that G dC'j , 

J a J « J a 

lies within D2.266 (see Figure |3)). For our choice of the curve 7, the return map of the set 

la J a J a 

is given by the 2-nd and the 3-d iterates of fa on L^, and -R°- , respectively. 

Ja Ja 

6. Hyperbolic properties of cylinder renormalization 

The expanding direction of TZcyi- It is not difficult to see that the operator TZcyi 
possesses an expanding direction at / (cf. |Ya3j ) : 

Proof of Theorem \2.4[ Part (III). Let v{z) be a vector field in Cu, 

v{z) = v'{0)z + o{z). 

Denote 7^ the quantity 

For a smooth family 
we have 

ft{z) = a';{z){f'{0)z + o{z)), where at(0) = 1 + ^7, + o{t). 
The g^+i-st iterate 

/r-(z) = ia^,iz)r--^iifiO)y--^z + oiz)). 

In the neighborhood of the renormalized vector field Cv is obtained by applying a uni- 
formizing coordinate 

Hence, 

«f (0) = [(ar(0))''-+f , 



m 

h{z) = f{z)+tv{z) + o{t), 
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SO 

Icv = ^Iv, where A = I3qm+i > 1- 

Hence the spectral radius 

Rs^iCv) > 1, 

and since every non-zero element of the spectrum of a compact operator is an eigenvalue, 
the claim follows. □ 

Numerical verification of hyperbolicity of TZcyi- It is natural to conjecture: 

Conjecture 6.1. There exists an open neighborhood U C Cu containing f such that TZcyi 
is a strong contraction in 

W = {f eU\ /'(O) = e'"'^}. 

Thus,W = W,l^if). 

To verify this conjecture numerically, we have to justify using a finite-dimensional ap- 
proximation to C to test for contraction. For this we rely on a numerical observation 
discussed in the previous section: 

£ : ^ Lj,, with p = 2.266, and p' = 3. 

This implies, that the finite-dimensional approximations of £ obtained by truncating all 
the powers higher that z'^ will converge geometrically fast in N. 
Set hj to be the coordinate vectors hj{z) — / p' , so that 

= sup ||>C/ij||p. 

Since a perturbation / -|- ehj docs not lie in L^, we perturb along a different set of vectors: 
(6.1) ^J = lAr^ g,{z) = z^ - ^z^+\ j>l. 



\9j\\p J + 1 



which form a basis in L^. 



Numerically, to estimate the spectral radius 



we can fix a large enough N and a small e (we have used the value e = 0.01), compute for 
each Cj, 2 < j < N, the finite difference 

1 

-{T^cylifa + eej) — TZcylfa), 

truncate past the A'"-th power - and expand over the basis vectors Cj to obtain an {N — 
1) X (A^ — 1) matrix Ajy- Below we present the approximate expression for Aq (the numbers 
have been truncated to the fifth decimal). 
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0.45879- 
0.97624i 


0.68789- 
0.46254i 


0.11338- 
0.09738i 


0.13041+ 
0.07490i 


0.15824+ 
0.11616i 


-0.13666+ 
1.72834i 


-0.64474+ 
0.54306i 


0.33937- 
0.50837^ 


0.14710- 
0.13849?: 


-0.21006+ 
0.02552i 


-0.90634- 
1.37322i 


0.27155- 
0.38270i 


-0.05765+ 
1.09081i 


-0.22948+ 
0.14700i 


0.18338- 
0.39078i 


1.23970+ 
0.20634i 


-0.08549+ 
0.23219i 


-0.68227- 
0.81153i 


0.12817- 

0.14685i 


0.10981+ 

0.47861i 


-0.63443+ 
0.58168i 


-0.02489- 
0.18893i 


0.80205+ 
0.00357i 


-0.04014+ 
0.12864i 


-0.34685- 
0.19936i 



Table 1. Matrix Aq. 

This matrix has the spectral radius 

i?Sp(A) -0.53. 

Estimating the spectral radius. We now proceed to produce a justification for the 
above numerical experiment. We will equip L^, viewed as a vector space, with a new 
Zi-norm 

oo 

(6.2) l/|p = El/^l' 

k=l 

where are the cocfiicicnts in the expansion of / in the basis {ej}: f = Yl'k'=i fk^k] and 
denote the new Banach space by L^. The projection P<jv on span^<^<^{ej} will be defined 
by setting 

N 

(6.3) J]/,.e,.. 

We will also abbreviate I — P<Ar as P>Ar. 

We would like to emphasize that P<Ar/ € Lp whenever / e L^, and therefore the operator 

(6.4) A = F<nCF<n 

serves as a finitc-dimcnsional approximation to the action of >C on L^. We will now make 
the latter statement more precise. 
To this end, notice, that 

(6.5) C = A + CF>N + FyNCF<N = A + n. 

The following Lemma demonstrates how one can obtain an upper bound on the spectral 
radius of the differential C at the fixed point in terms of the norm of a power of the 
finite-rank operator A and the magnitude of the norm of H. 
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Lemma 6.2. Let C = A + 7i he a hounded operator on some Banach space, such that 
^A^W < 7 < 1 for some k > 1 and, \\Ti\\ < 5 < 1. Then, the spectral radius Rsp{C) 
satisfies 

(6.6) Rs,{C)<^''\l + C5hY'^ 

for some (explicit) constant C. 

Proof. The claim follows from the spectral radius formula. First, 



Rsp{C) = lim^^ooll-^"!!" = lii^n- 



< lim„_ 



lim^ 



fciii- 



The norm ||>C'^|| is finite, and therefore 



lim„^oo||>^^ 



Then 



Rsp{^) < linin- 



< lim^ 



lim^^ooll-^ 11 < lim^ 



Let nCk denotes the binomial coefficients, and let C = Yli=i kCi\\A'' ""^ then 



Rsp{jO-) < lim„_ 



WA'^r + J2nCi\\AT-'{c5y 



i=l 



< hm^^oo [7"(1 + ((C(5/7 + 1)'^ - 1)]^ 



□ 



It is left now to bound the difference of JC from A. First, we state the following Cauchy- 
type estimate, whose straightforward proof will be left to the reader: 

Proposition 6.3. Assume that an operator TZ^yi is analytic in an open hall Bj.{f ) C L],. 
Let e < 1, and h Ehp be such that \h\p < r. Then 

(6.7) \T^cyi{f + ^h) - Ucyif - eCh\p < sup \TZcyi{f + sh) - Tlayiflp. 

J- - e |s|<i 

Note that 



mTfw\p < sup \£ej\p. 



This, together with the preceding Proposition and the compactness property of renormal- 
ization, immediately implies that can be bound by a finite difference. Specifically, 
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for all j > N: 

e"^ sup \TZcyi{f + ^h) - f\p + sup \TZcyi{f + h) - 



hlp<l l'llp<l 

e 



< e-^ sup |P<Jv[7^cJ,/(/ + e/^) -/]!, + -— sup \f<N[n,yi{f + h)-f]\p 

/i|p<l |h|p<l 

' \ N+1 

4 sup |P>^[7^,,K/ + e/^)-/]|p' + ^ sup |P>;v[7^c,K/ + /^)-/]| 



|(i|p<i l'i|p<i 



(6.8) = Ci, 

Similarly, for all 2 < j < N: 

\^>NCe,\p < e-^ sup \FyN[n,yi{f + eh)-f]\p + -^ sup |P>^[7^e^,^(/ + Z^) - /II 



heTjtW ' i — 6 heT-w 

\p<l \h\p<l 

e 





N+l 








I sup ^ 




/ heTj:W 




lh\p<l 


C2. 





\h\p<l 

(6.9) 

We would like to emphasize that these bounds use the fact that £ is a compact operator 
in an essential way. Bounds (|6.8p and (|6.9|) can be used to estimate \{C — ^)|T^.i^|p. In 
particular, according to equation ()6.5p . 

|(£ - A)\Tfw\p < max ||£P>;v|p, \F^nC\t^w^<n\p] < max{Ci, C^} . 

This expression provides a bound on 6 in 

We have chosen N = 14, and experimentally bounded Ci and C2 by testing on vectors 
h = ei5, and h = 62 which empirically maximize the respective suprema. Choosing A; = 80 
in ()6.6p . we have the following values for the constants that enter estimate ()6.6|) : 

7 < 2.07 X 10"l^ (5 < 0.24, C < 8.4 x 10-^ 

therefore, according to 

i?Sp(>^:|T,-w^) < 0.85. 

A better bound can be obtained if one computes all relevant constants for a larger value 
of N, which requires more computer time. It is plausible that the spectral radius is close 
to 0.58, since we have observed that as increases the largest eigenvalue of the operator 

converges to 

A = 0.15 + 20.56. 
This eigenvalue has been truncated to 2 decimal places. 
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As a final comment, note that the following simple observation implies that perturbations 
in the directions of the vectors hj can also be used for estimating the spectral radius of 

Proposition 6.4. We have 

Spec{C\-Bi) = Spec{C\-Li). 

To see this, note that the only difference between the spectra is that (contained in both 
spectra) is an eigenvalue of the operator £|bi corresponding to hnear rescalings. We leave 
the straightforward details to the reader. 



7. Appendix 

The objective of this Appendix will be to describe how the Cauchy (j3.8j) and Hilbert 
()3.7j) transforms can be computed numerically. 

The Constructive Measurable Riemann Theorem 13. II deals with Lp functions which gen- 
erally do not need to be different iable. Therefore, one has to chose an appropriate repre- 
sentation of the Lp functions that enter the Theorem 13. H possibly, as a collection of values 
on a grid, or as a Fourier series with the radially dependent coefficients. The latter choice 
has been made, for instance, in |Dalj . |Gaij . |GKj . and will be also adopted in the present 
paper. 

Represent h and P[h] in (|3.8|) as: 

oo 

(7.1) h{re'') = Yl hk{r)e''', 

k=—oo 
oo 

(7.2) P[h]{re^') = ^ Pk{r)e^'', 

k=—oo 

where the coefficients of the P-transform are given by 

(7.3) pk{r) = ^ r e-^^'P[h]{re'') dO. 

A classical theorem of analysis (cf jAhlj . |Mar j ) states that Cauchy transform of an Lp- 
function, p > 2, is well-defined and is Holder continuous with exponent 1 — 2/p. In |Dalj 
and |GKj this fact has been used to show that the Fourier coefficients of Cauchy transform 
are given by the following equations 

2^ Q hk+iip)dp, k<0, 



(7.4) pk{r) 



-2^ (^^^ hk+i{p)dp, k>0. 



To obtain similar formulae for Hilbert transform, assume that h is a Holder continuous 
function compactly supported in an open disk around zero of radius R, B{0, R) C C. The 
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Hilbert transform of such function is known to exist as a Cauchy principal value (cf |Ahlj . 
Jarj ) ■ As with Cauchy transform, represent this transform as a Fourier series: 



(7.5) T[h]{re'') = V Cfc(r)e^'=^ Cfc(r) = — / e-'''T[h]{re'') dO 



In |Da2j , |Da3j and |GKj the authors arrive at the following expressions for these coeffi- 
cients: 

(7.6) co(0) = -21im ! ^^^^dp, and Cfc(O) = 0, whenever ^ 0, 

(7.7) Cfc(r) = Ak -—jhk+2ip)dp + Bk / -—^hk+2{p)dp + hk+2{r), 

Jo P Jr P 

where 

^ ^ " \2(A; + 1) ,A; < 0, \ 0, A; < 0. 

We would like to mention that the fact that Hilbert transform is a singular integral op- 
erator makes a rigorous justification of the formulas ()7.6j) - ()7.7jl significantly more involved 
than that of (fT^ . 

Formulas ()7.4|1 and ()7.(i|l - ()7.7|l can be used to construct an efficient algorithm for solving 
a Beltrami equation. Given values of h, for instance, on a circular N x M grid that 
contains the compact support of h, one can use a fast Fourier transform (FFT, cf |NRj ) 
to find the values of the coefficients hk at the radii Ti, 1 < i < M. Next, one can use 
these values to construct a piecewise constant, a piecewise hnear or a spline approximation 
of the functions hk (the choice of approximation, of course, depends on the known or 
expected smoothness of hk). This allows one to compute integrals in ()7.4|1 and in ()7.6|) - 
(j7.7|) . Armed with these implementations of Hilbert and Cauchy transforms, one can try 
to solve the Beltrami equation (j3.3|) . first by running iterations (j3.1(J|) for some time, and, 
finally, applying (j3.1H) . It is convenient to use the point- wise multiplication of grid values 
of h and p + 1 inside Hilbert transform in (j3.1Up . rather than the multiplication of their 
Fourier series: the order of the computational complexity of the point-wise multiplication 
is 0{NM), as opposed to 0{NM'^) for the series. The transition from the representation 
of /i as a Fourier series to point values at each iteration step can be performed with the 
help of the FFT. This way, the computational complexity of one iteration step becomes 
0(iVM log2 M). 
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